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Nonexponential relaxation of magnetization at resonant 
tunneling points of nanoscale molecular magnets is inter- 
preted to be an effect of fluctuating random field around the 
applied field. We demonstrate such relaxation in Langevin 
equation analysis and clarify how the initial relaxation 
(square-root time) changes to the exponential decay. The 
scaling properties of the relaxation are also discussed. 

PACS number: 75.40.Gb,76.20.+q 

As to the relaxation of metastable magnetization of 
uniaxial nanoscale molecular magnets such as Mni2 and 
Fes, the resonant tunneling phenomena have been paid 
attention and various interesting properties of the phe- 
nomena have been reported [^-|9| . Properties of the phe- 
nomena have been also investigated in the theoretical 
side pO|-pO|. We have studied such relaxation from a 
view point of nonadiabatic transition among their dis- 
crete energy structure due to a finite number of degrees 
of freedom. There we studied the relation between the 
relaxation rate and the sweeping velocity. Actually such 
relation has been utilized to estimate the energy gap in 
a recent experiment [§J^|. Because the tunneling gap is 
so small, thermal environments have strong influence on 
the tunneling, which has been studied extensively. 

In the case of static field, there have been many exper- 
iments reporting that nonexponential relaxation occurs 
at the resonant point. That is, the initial relaxation can 
not be fitted by usual exponential function but is fitted 
well by a stretched exponential function with the expo- 
nent near 0.5 or a square-root function [Q(2^-|2^]. If the 
field would be precisely fixed at the resonant point, we 
should see the coherent tunneling. Thus the observa- 
tions indicate that the decoherent effect plays an impor- 
tant role. This nonexponential relaxation has been inter- 
preted as a phenomenon due to the fact that the region 
of the field of the resonant tunneling is very narrow and 
small fluctuation can detune the resonance condition by 
Prokof 'ev and Stamp [Q . They considered the distribu- 
tion of the internal field which mainly consists of dipolar 
field from other molecules, and investigated its time evo- 
lution. Combining the evolution of the distribution and 
the relation between the steady state distribution of the 
dipolar field and the magnetization, they found a square- 
root time initial relaxation and an exponential relaxation 
in the late stage, which explains the overall dependence 
of the experiments. 



In the detailed observation on the distribution of in- 
ternal field P{£) in Fe§ M, the square-root dependence 
is found even if the initial magnetization is zero, and 
change of the distribution appears only locally and it is 
not associated with the reforming total distribution of 
P(£). These features are not compatible with the mean- 
field type analysis of P(£), and thus it seems necessary 
to consider more general mechanism of the square-root 
time initial relaxation. Actually it has been pointed out 
that the fast fluctuation of hyperfine fluctuations are 
very important in the local field dynamics. In this Letter, 
we propose an alternate explanation of this phenomenon 
as a general phenomenon at narrow resonant points with 
fluctuating field using a Langevin equation approach, i.e., 
using the Ornstein-Uhlenbeck process f25|j . 

Here we consider a two-level system for the simplicity 
and its Hamiltonian is given by 

AP 

H = h cxt (t)a z - —a x , (1) 

where a x and u z are the Pauli matrices and h ext (t) is a 
time-dependent external field, and 4p denotes the trans- 
verse field which represents the quantum fluctuation of 
the z component of the spin. We solve the Schrodinger 
equation of this system with time dependent field h cxt (t) 
by applying the time evolution function 

\t + At) =exp(-iU(t))\t), (2) 

where the h ex t(t) is prepared by a Langevin equation. 

Corresponding to the molecules of Mni2 or Fes, we 
consider that the resonant point is very narrow and is 
regarded as a point. The external field consists of the 
static part ho and a fluctuating part h(t): 

h eKt (t) = h + h(t). (3) 

The fluctuating part is caused by independent changes of 
many magnetizations around the site. Thus we consider 
that this part is an assembly of independent fluctuation 

h(t) = h(0) + f y2 5h i( s ) ds > ( 4 ) 
Jo j 

where Shj (s) is the change of a field from the j-th site (hj) 
at time s. Here we assume that Shj(s) is independent of 
time and position, and we regard it as a white gaussian 
noise rj(s) 



1 



(V(t))=0, (v(t)v(s))=2DS(t-s). 



(5) 



Here the assumption of the gauss distribution is not es- 
sential but we assume it for the convenience of analysis. 
Then, the relation^]) is written as 



h(t) = h{0) + / i](s)ds = W(t), 
Jo 

where W is the Wiener process 

(W(t)} = 0, (W(t)W(s)) = 2DMin(i, s). 
The distribution of W(t) is given by 
1 



P{W) 



exp 



W 2 
4Dt 



(6) 



(7) 



(8) 



We show an example of the Wiener process in Fig. 1 (a 
dashed line). 

The nonadiabatic transition only occurs when the field 
crosses the resonant point. Hereafter we set the reso- 
nant point at h = and also the static field at ho = 0. 
The probability p of the adiabatic transition of the state, 
i.e., from the ground state to the ground state (or from 
the excited state to the excited state), is given by the 
Landau-Zener-Stiickelberg (LZS) formula [|26|- 



p = 1 — exp — 



(AE)< 



2\M-M'\v 



(9) 



where AE is the energy gap at the crossing point (tun- 
nel gap), M and M' are the magnetizations of crossing 
states, and v is the speed of the field at the resonant 
point. The transition probability from the initial distri- 
bution (j>i,P2) to the scattered distribution (p'l,^) de- 
pends on the velocity v of the field at the time of crossing 
and also on the phase factor due to the free motions out- 
side of the crossing region p9|| . Here we consider many 
sample of {h cxt (t)} and the ensemble average over the 
distribution of v and the phase factor. 

The change of magnetization (M — * M') at this cross- 
ing is given by M = M {pi - p-i) -» M' = M^(p[ - p' 2 ), 
where Mq and Mq are the magnetization of the each 
ground state. In Fig. 1, we show a dynamics of mag- 
netization M(t) for the shown process of h eX f 

From the dependence of (|9|), the average probability of 
the adiabatic transition, which changes the magnetiza- 
tion at a crossing, is estimated as 



Pa. 



l 2 M -M' 



a - 



(A£) 2 



(10) 



where we assumed that the value of p is small and that 
the average velocity of the field is proportional to the 
strength of the random field \[~D. The change of magne- 
tization at a crossing is given by M' = (1 — 2p av )M. 




i» 1 \ 



2000 



FIG. 1. A sample of the Wiener process (dashed line) and 
the change of the magnetization for this process of /i ex t ■ 

If we assume that the sequence of the crossings give in- 
dependent contributions, the change of magnetization for 
arbitrary initial value M (0) is given by 



M(t)=M(0)(l-2p av ) JV -«, 



(11) 



where N(t) ac is the accumulated number of the crossings 
by the time of t. 

The number of crossings per unit time, N(t), is esti- 
mated from the probability for \W\ < VD, where \f~D is 
the jump range of the field. Thus N(t) at time t is given 
by 



N(t) oc —== x 4~D oc 



(12) 



Thus the accumulated number of crossings by the time t 
is given by 



iVac(«)= / N(t)dt = cVt, 



(13) 



where c is a constant. Thus the total change of the mag- 
netization is 



M = M (l - 2a Q 



exp 



(14) 



where a = 2a^c. This constant is uniquely determined 
by the nature of random process which will be deter- 
mined later. Thus we conclude that the magnetization 
shows a stretched exponential decay with the exponent 
1/2, which shows the square-root time initial relaxation. 
In Fig. 2, we show the averaged magnetization (M(t)) 
over the 10,000 samples of h cxt . In the inset M(t) is 
plotted in the coordinate U/t, log M(t)), where we con- 
firm the dependence of (|l4|). This dependence gives the 
initial square-root time dependence 
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M = M (l 



(15) 
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FIG. 2. The averaged magnetization (M (t)) over the 
10,000 samples of /i e xt- The inset shows a plot (log M(t), t 1 ' 2 ), 
where the dotted line is a straight guide line for eye. 

In the real situation the field h(t) may not deviate 
without limit and there are some restoring mechanisms. 
Next we study the mechanism of transition from the ini- 
tial relaxation (|l4|) to the exponential relaxation. Taking 
into account the restoring effect, we adopt the Ornstein- 
Uhlenbeck process for the evolution of h(t): 



dh 

~dt 



= - 1 h(t)+T)(t). 



The process (g) corresponds to the case 7 = 0. 
distribution of h(t) for this process is given by 



P(h) 



1 



V2 



: exp 



where 



D 



2a 2 



cr = — [1 - exp(-27t)] 
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(16) 



The 



(17) 



(18) 



For t <C 7 _1 , this reproduces ffi). On the other hand 
it has a stationary distribution a 2 = D/^f at large time 
t ^> 7 _1 . At this late stage there is a constant probability 
that the h(t) stays around 0, which causes a constant rate 
relaxation, i.e., the exponential relaxation. Therefore we 
expect 



M = Mi exp 




(19) 



for a long time t ^> (27)" 1 , while we have the relation 
( |l4| ) for a short time t <C (27)~ 1 . Here (3 is a constant 
independent of T, 7, and D. Mq and Mi are constants 
corresponding to a kind of initial magnetization of each 
process, which does not coincide with the initial magne- 
tization. 



In Fig. 3(a) we show a time dependence of the mag- 
netization (M(t)) averaged over 10,000 samples. In Figs. 
3(b) and (c), we plot M{t) in {s/i,\ogM{t)) and in 
(t,\ogM(t)), respectively. There we find the crossover 
from the square-root time relaxation to the exponential 
relaxation around t ~ (27) -1 . 
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FIG. 3. Parameter dependence of (M(t)) for /i ex t given by 
@ for D = 0.001: (a) the linear plot (t,(M(t))), (b) plot 
(yt, log(M(t))) for the initial relaxation (14]), and (c) plot 
(i,log(M(t)}) for the late relaxation (|l|). 

In a very early time t < 5 we find a dead time where 
the magnetization does not change. This is considered to 
be due to the fact that the noise starting at stays near 
the resonant region where the system essentially evolves 
coherently, i.e. M(t) ~ cos(AEt). This dead time be- 
comes relatively short when AE and 7 becomes small. 

Studying (M(t)) for various parameters (AE : ^,D) 
we confirm that scaling relations of M(t) on the pa- 
rameters indicated by ( fl4"| ) and jl9|). In Figs. 4(a) 
and (b), we plot the data (AE — 0.01, 7 = 
0.001,0.002, and 0.003, and D = 0.001,0.002 and 0.003) 
in the coordinates: (V~t,\ogM(t)/((AE) 2 fy/D)) and 
(((AE) 2 y^jD)t,\ogM(t)), respectively. We could plot 
the data in (((AE) 2 /y/D)Vi, log M(t)) instead of Fig. 
4(a). However the scaling time region [the dead time ~ 
0(1) < t < (27) -1 ] of each parameter set shown in dif- 
ferent region of ((AE) 2 /^D)Vt and it looks mess. Thus 
we plot data in the way of Fig. 4(a). There the lines 
shows some distribution due to the distribution of Mq 
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and Mi, and also the effect of the dead time. However, 
we find that slopes of them are almost the same to each 
other and we estimate the constants as 



q - 3.5 ±0.1, and O ~ 2.3 ±0.1. 



(20) 



When we change AE, the value of Mq and Mi and the 
dead time change, but we find that the values of a and 
(3 are consistent although the scattering of date is larger. 
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FIG. 4. Scaling plots (a)(v^, log M {t) / {AE 2 / ^D)) , and 
(b)(((AE) 2 yfi/D)t, log M(t)) 

Finally we consider the case where the property of the 
noise depends on M(t). Generally the change of the av- 
erage magnetization causes the change of distribution of 
the field /i ex t as discussed by Prokof'ev and Stamp p4| . 
In our analysis this change should be taken into account 
as a slow change of ho, which may lead the same mech- 
anism discussed in Jlif , which will be discussed in the 
future. Effect of change of the mean field to the LZS 
transition has been also discussed as a feedback effect 
which causes a large change of the transition probabil- 
ity 1 30 1 . This effect also causes important modification 
when the field is swept where the effective sweeping rate 
is modified. In the present case of fixed field, we assume 
that the time scale of the fluctuation field is much smaller 
than the change of the average magnetization. 
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